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fvj ' Using the Schwinger-Keldysh technique, we derive the transport equations for a system of quantum 

, scalar fields. We first discuss the general structure of the equations and then their collision terms. 

M ■ Taking into account up to three-loop diagrams in (j)^ model and up to four-loop diagrams in 0'* 

model, we obtain transport equations which include the contributions of multi-particle collisions 

and particle production processes, in addition to mean-field effects and binary interactions. 
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I. INTRODUCTION 



> 

a^ 

^^ ' Transport theory is a very convenient tool to study many-body nonequilibrium systems, both relativistic and 

>P , nonrelativistic. The kinetic equations, which play a central role in the transport approach, usually assume that 

dissipation processes are governed by binary collisions. However, when the system of interest is very dense one 

expects that multi-particle interactions will play a significant role. Such interactions are known to the control spectra of 

fluctuations and transport properties of dense gases and liquids |l| . Furthermore, in relativistic systems a characteristic 

i-Q particle's kinetic energy is usually comparable to the particle's mass, and processes leading to particle production 

^^' become important for a system's dynamics. 

O i' It is therefore expected that multi-particle interactions and production processes will play an important role in the 

D , dynamics of a relativistic quark-gluon plasma of high energy density. The importance of gluon multiplication in the 

i-C ' process of the plasma's equilibration has been repeatedly stressed, see e.g. j^., 3, ^, ;&, ^sJIJ- The scattering of three 

^ , gluons into three gluons has been studied in the context of thermalization in |^ . It has also been shown within the 

; '"! ' scalar field theory that the transport coefficient of bulk viscosity, as given by the Kubo formula, strongly depends on 

rS \ particle number changing processes 9, 10]. We conclude that a complete description of the quark gluon plasma based 

jrt ' on transport theory requires the derivation of the relativistic transport equation which includes the multi-particle 

■ " " ' collisions and particle production processes, in addition to the mean-field effects and binary interactions. The form 

of such an equation has been postulated by some authors, see e.g. 6, Sl,y,[l3, but a systematic derivation from first 

principles is lacking even in the simplest scalar field theory. 

The Schwinger-Keldysh [l^, u2^ formulation of quantum field theory provides a framework for the derivation of the 
transport equation. Kadanoff and Baym [13 developed a technique for nonrelativistic quantum systems which has 
been further generalized to relativistic systems [IlllllllllllllllllHllllilliU^. We mention here only the 
papers which go beyond the mean field or Vlasov approximation and provide a more or less systematic analysis of the 
(binary) collision term. 

In this paper, we consider self-interacting scalar fields with cubic and quartic interaction terms and study their 
transport equations beyond the binary collision approximation. There are basically two separate parts of our study. 
The first part is the derivation of the general structure of the transport equations. Such an analysis has been published 
previously by one of us [lalSfllsS and it is presented here only to provide a framework for the second part of the paper. 
The main steps of the derivation are the following. We define the contour Green's function with the time arguments 
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on the contour in the complex time plane. This function is a key element of the Schwinger-Keldysh approach. After 
discussing its properties and relevance for nonequilibrium systems, we write down the exact equations of motion i.e. 
the Dyson-Schwinger equations. We perform a Wigner transformation and do a gradient expansion by assuming that 
the system has macroscopic quasi-homogeneity. The resulting pair of Dyson-Schwinger equations is converted into the 
transport and mass-shell equations which are satisfied by the Wigner function. We define the distribution function 
of usual probabilistic interpretation, and we find the transport equations satisfied by this function. The transport 
equation is derived to lowest order in the gradient expansion, which corresponds to the Markovian limit in which 
memory terms are neglected. 

In the second part of our study, we perform a perturbative analysis of the self-energies which enter the transport 
equation. We show that the Vlasov or mean-field terms are dominated by the lowest order tadpole diagrams, while 
binary collisions emerge from two-loop contributions. Multi-particle interactions and production processes appear at 
three-loop level in the cj)'^ model and at four-loops in the cj)'^ model. There are obviously higher loop contributions to 
both mean field and binary collision terms but these higher loop contributions will not change these effects qualitatively, 
and for this reason are not studied here. Instead we obtain leading order contributions to the collision term that result 
from 2 <-!■ 3 processes in (j>^ model, and 2 <-^ 4 and 3 ^^ 3 processes in (j)* model. To extract the processes represented 
by the three- and four-loop graphs, we work with the Keldysh representation, using a technique developed previously 
by one of us pll27t. 

Multi-particle and production processes have already been discussed in the Kcldysh-Schwinger approach to transport 
theory but, to our best knowledge, only for nonrelativistic systems ^8, 29, 30J. In this case the analysis is quite 
different. In the nonrelativistic limit, interactions are instantaneous and internal lines representing the potential 
interaction cannot be cut. Consequently, the extraction of the scattering matrix elements in the small coupling limit 
is considerably simpler than in the relativistic case. However, it should be noted that much progress has been achieved 
in resuming multi-loop diagrams in the nonrelativistic approach [29j. 

The problem of off-mass-shell transport has recently attracted a lot attention in the literature |3lL 133 . l3a. l34L l35ll3q | . 
The collision term for off-mass-shell transport contains important contributions from one- and two-loop self-energy 
diagrams which correspond to 1 ^^ 2 and 1 ^^ 3 processes. These processes are kinematically forbidden for on-mass- 
shell particles, and for this reason are not studied here. 

Throughout this article we use natural units where h = c = 1. The signature of the metric tensor is (+,—,—,—). 

II. PRELIMINARIES 

We consider a system of real scalar fields with Lagrangian density of the form 

C{x) = ^d^'^{x)d^^{x) - \m'<i?{x) - ^rix) , (1) 

where n equals 3 or 4. The renormalization counterterms are omitted here. The field satisfies the equation of motion 
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[d^ + m^]^{j:)^-^-^^--\a:). (2) 



The energy-momentum tensor is defined as 

T^'-^ix) = d^'(l}{x)d''(t)ix)-g'"'Cix) . 
Subtracting the total derivative 

^d^'d''cj)^{x)-g^^'^^d''d„cj)\x), 

and considering the case of free fields, we obtain an expression for the energy-momentum tensor which has a form 
that is convenient for our purposes 

Tr{x) = -^^ix)TT Hx). (3) 

III. GREEN'S FUNCTIONS 

We define the contour Green's function as 

^A{x,y)''^'{m^my)), (4) 



where (...) denotes the ensemble average at time to (usually identified with — oo) and T is the time ordering operation 
along the directed contour shown in Fig.^ The parameter imax is shifted to +(X) in calculations. The time arguments 
are complex with an infinitesimal positive or negative imaginary part which locates them on the upper or lower branch 
of the contour. The ordering operation is defined as 

f(t){x)(t){y) = e{xo,yo)(f>{x)(f>{y) + e{yo,XQ)(j){y)(p{x) , (5) 

where 6(xo, j/o) equals 1 if xq succeeds yo on the contour, and if a;o precedes j/q. 

1-max 



FIG. 1: The contour in the complex time plane for an evaluation of the operator expectation values. 

If the field develops a finite expectation value, as is the case when the symmetry is spontaneously broken, the 
contribution (0(a;))(0(2/)) is subtracted from the right-hand-side of the equation defining the Green's function and 
one looks at field fiuctuations around the expectation values (see e.g. [ial23|)- Since (i^(a;)) is expected to vanish in 
the models defined by the Lagrangians ^ we neglect this contribution in the definition of the Green's function. 

The contour propagator involves four functions, which can be thought of as corresponding to propagation along the 
top branch of the contour, from the top branch to the bottom branch, from the bottom branch to the top branch, 
and along the bottom branch. Below, we define four Green's functions with real time arguments, and discuss their 
relationship to the contour Green's function. We define 

tA>{x,y) 1^' (0(x)0(y)), (6) 

tA<{x,y) "^^ {cb{y)cf>{^)) , 

zA^(x,y) =' (T=0(x)0(2/)) , 

tA'^{x,y) 1i' (T>(.T)0(y)), 
where T'^{T°') indicates chronological (anti-chronological) time ordering 

T''(j){x)(j){y) = e(a;o - yo)(p{x)(p{y) + Q{yo - xa)(t){y)(j){x) , 

rXx)0(2/) 1^' 9(2/0 - xo)0(a;)0(y) + e(xo - yo)c^{y)<l>{x) . 

These functions are related to the contour Green's functions in the following manner: 

A!^ [x , y) = A{x , y) for a;o , yo on the upper branch, (7) 

A°(a;, y) = A(a;, y) for xq , yo on the lower branch, 

A^(x, y) = A(a;, y) for xq on the upper branch and 

yo on the lower one, 

A^(x, y) = A(a;, y) for a;o on the lower branch and 

yo on the upper one. 

The top branch of the contour is usually called the '1' branch and the bottom branch is called the '2' branch. Using 
this notation, the functions A"^, A°, A< and A> are often written as 

A^ =. An , A'^ = A22 , (8) 

A< = A12 , A> = A21 . 



The following identities can be obtained directly from the definitions given above: 

A"(a;,y) = Qix„ - y„)A^{x,y) + Q{yo - Xo)A^{x,y) , (9) 

(iA°(x,2/))t = iA%x,y) = iA^iy,x), 

(iA^(x,2/))t = iA^{x,y)^iA^{y,x), 

where the f denotes hermitian conjugation which involves an interchange of the arguments of the Green's function. 
Using Eqs. ® and @, it is easy to show that the four components of the contour Green's function are not independent 
but satisfy the relation: 

A%x,y) + A'^{x,y) - A<{x,y) - A>{x,y) = (10) 

In some situations, it is useful to work with retarded (+), advanced (— ) and symmetric Green's functions. We 
make the usual definitions for these propagators: 

iA+{x,y) 1i' e{xo-yo)mx),cl,{y)])^iQ{xo-yo){A>{x,y)-A<{x,y)) , (11) 

iA-{x,y) = ~e{yo-xo)m^),^{y)])^teiy,^xo){A<{x,y)-A>{x,y)) , (12) 

iA^y'^ix,y) 1^*' {{cj,ix),(j,iy)}) ^ i{A>ix,y) + A<{x,y)) , (13) 

where the curly brackets indicate an anticonrmutator. 

It is easy to show that there is a simple relationship between the original set of propagators (A'^, A°, A< and 
A>) and the retarded, advanced and symmetric propagators (A+, A~ and A'^y™). In fact, since the four original 
propagators satisfy the constraint (|10|l and thus only three of them are independent, transferring between the two 
representations is equivalent to a change of basis. Henceforth, we will refer to the set of propagators A+, A~ and 
^sym g^g ^]-^g Keldysh representation of the contour ordered propagator, and the set A'^, A°, A< and A> as the 1/2 
representation. Using the definitions © and the relations (O, it is straightforward to connect the Keldysh to the 1/2 
basis 

A^{x,y) = A'{x,y)-A^{x,y), (14) 

A^(x,y) = i(A^y-(x,2;)±A+(x,y)TA-(x,2/)) , 

A+(x,y) - A-{x,y) ^ A>{x,y) ^ A<{x,y) . 

The relations H9I10I14|I hold for real and complex fields. For pure real fields, which are studied here, there are extra 
relations among Green's functions. In particular, we will make use of the idendity 

A>ix,y)^A<iy,x). (15) 

We discuss briefly the physical interpretation of the Green's functions we have defined above. The functions A° and 
A describe the propagation of a disturbance in which a single particle or antiparticle is added to the many-particle 
system at space-time point y and then is removed from it at a space-time point x. The function A''{x,y) describes 
a particle disturbance propagating forward in time, and an antiparticle disturbance propagating backward in time. 
The meaning of A°(x,y) is analogous but particles are propagated backward in time and antiparticles forward. In 
the zero density limit A'^{x, y) coincides with the Feynman propagator. In the case of the retarded (advanced) Green 
functions, both particles and antiparticles evolve forward (backward) in time. 

The physical meaning of functions A^(a::, y) is more transparent when one considers their Wigner transforms which 
are given by 

A^{X,p) l^f J d% e^P«A^(X + ^u,X~ ^u) . (16) 

It is easy to show that the free-field energy-momentum tensor (O averaged over the ensemble can be expressed as 

(To'-^iX)) = J ^p>'^*A^(X,p) . (17) 



One recognizes the standard form of the energy-momentum tensor in kinetic theory with the function zA^(X,p) 
playing the role of the density of particles with four-momentum p at the space-time point X. This observation leads 
us to consider i A^ {X, p) as the quantum analog of the classical distribution function. The function iA^ {X, p) is 
hermitian but it is not positive definite, and thus the probabilistic interpretation is only approximately valid. One 
should also note that, in contrast to the classical distribution function, the function iA^(X,p) can be nonzero for the 
off-mass-shell four-momenta. As will be explicitly shown in Sec. IIXI one can extract the usual distribution function 
from iA^(X,p) under specific conditions. 

IV. EXACT EQUATIONS OF MOTION 

Since the Green's function j A^ [X, p) will be used to obtain the distribution function, we derive an equation which 
governs the evolution of zA^(Ar, p). We start with the Dyson-Schwinger equation for the contour Green's function 
which has the form 

A{x, v) - Ao(x, y) - f d^z f d'^z' Ao{x, z)n{z, z')Aiz', y) , (18) 

Jc Jc 

where Ao(a;,y) is the free Green's function. We note that the integrals over zq and Zq in the last term of equation 
(|18|l are performed along the contour. 

The Dyson-Schwinger equation (|18l) can be rewritten as two equations: 

[dl + m^] A(x, y) = -<5(4) (x, y) + f d'x'n{x, x')A{x\ y) , (19) 

Jc 

[dl + m"] A(x, y) = -5(4) (^^ y^^ f rf4^/^(^^ ^/)n(x', y) , (20) 



c 



where the function 6'^^\x, y) is defined on the contour as 



{(5*-4' [x — y) for xq , j/o from the upper branch, 
for xq , j/o from the different branches, 

—(5(4' {x — y) for xq , j/o from the lower branch. 

To obtain equations for A^ we split the self-energy into three parts as 

n[x,y) = ns{x)5^^\x,y) + n>[x,y)Q{x^,y^) +\l<{x,y)Q[ya,x^) , (21) 

where II^ corresponds to the tadpole contribution to the self-energy. With the help of the retarded and advanced 
Green's functions (|11I12|I and similar expressions for the retarded and advanced self-energies, we can rewrite equations 
(|19I20|I to obtain 

[dl + m''-ns{x)\A^{x,y) = j d^x'\n^{x,x')A- {x' ,y) + X]+ {x,x')A^{x' .y)^ (22) 

[92 + m2-n,(y)]A^(x,2;) - j d^x'[A'^{x,x')n- {x' ,y) + A+{x,x')n^{x' ,y)\ , (23) 

where all time integrations run from — oo to -|-cx). 

Similarly, the Dyson-Schwinger equations (|19I20|I provide the equations satisfied by the retarded and advanced 
functions A^(a;,j/) as 

[dl + m''-Us{x)\A^{x,y) = -S^^^x - y) + f d^x'U^{x,x')A^{x' ,y) , (24) 

[dl + m^-Ilsiy)]A^ix,y) = -5^^\x - y) + j d'x' A^{x,x')n^{x' ,y) . (25) 

We note that the equations for the retarded and for advanced functions are decoupled from each other - there is no 
mixing of the retarded, advanced and symmetric components of the propagators. 



V. APPROXIMATIONS 

It is a very difficult task to study systems wliicli are strongly inhomogcneous or strongly interacting - there are no 
general methods applicable to such systems. Thus, we assume that systems of interest are weakly inhomogcneous and 
weakly interacting. We discuss below the mathematical formulation of the corresponding approximations. 

1) The system is assumed to be weakly inhomogcneous in comparison with two different length scales. 

(a) We assume that the inhomogeneity length is large compared to the inverse of the characteristic momentum 
i.e. 

where F is either the propagator or the self-energy. This assumption allows one to perform the gradient 
expansion. 

(b) The inhomogeneity length of the system is assumed to be large compared to the inverse characteristic mass 
of free quasi-particles (or the Comptom wavelength) : 

\F{X,p)\ » l-\-JL—F{X,p)\ . (27) 



m 



dX^dXf' 



"-M 



This condition justifies the so-called quasi-particle approximation. When the bare fields are massless or the 
free mass is much smaller than the dynamically generated effective mass m*, the mass in Eq. (|27|l should 
be replaced by to* [23 . 

2) The assumption that the system is weakly interacting means that the coupling constant g in Eq. ^ is small 
and that all self-energies can be expanded perturbatively in g. 

The conditions of weak inhomogeneity and weak interaction will be used in five separate ways: 

(i) We will convert the equations (|22l I23|l into transport equations by performing Wigner transformations Hlti|) on 
all Green's functions and self-energies. Using Eq. (|26|l . we obtain the following set of translation rules: 

d^x'fix,x')gix',y) -^ f{X,p)g{X,p) + '-{f{X,p), g(X,p)} , (28) 

h{x)g{x,y) ^ h{X)g{X,p)~^-^^ d^ ' 

Hy)gix,y) ^ HX)giX,p)+^y-h^'-^, 

dl^f{x,y) -^ {-ip>^ + \d^)f{X,p) , 

d^fix.y) -^ {ip'^ + ^dnf{X,p), 

where we have introduced the Poisson-like bracket defined as 

dC{X, p) dD{X, p) dC{X, p) dD{X, p) 



[c{X,p), D{X,p)] 



dp^, 9X^ 9X^ dp^. 



The function h{x) is weakly dependent on x and we use the notation 9^ = n4— • This calculation is done in 
Sec. EU 

(ii) The condition ()27|l allows one to drop the terms containing d^ . These terms and those involving Imn+, which 
are discussed in (iv), are responsible for off-mass-shell contributions to the Green's functions A^. 

(iii) The assumption (|26|l is used together with the condition of weak interaction to show that the gradient terms in 
the right hand side of the transport equation (|35|l and mass-shell equation (|35|l can be dropped. This point is 
discussed in detail in Sec. IVII 



(iv) The weakly interacting approximation is used in Sec. |3 For 1/g :^ 1 and a finite bare mass, we have 

m^ - Us - Rcn+ > |Imn+| . 



(29) 



This assumption combined with the condition H27|l allows us to obtain an expression for the spectral function 
that is proportional to S{p^ — rn^). We note that both of equations H27(l and H29(l are commonly called 'the 
quasi-particle approximation'. This terminology is a consequence of the fact that there are two independent 
means by which the quasi-particle picture can be destroyed: strong inhomogeneity and strong interactions. 

(v) The weakly interacting condition is also used in Sec. IXIII to perform a perturbative expansion of the self- 
energies in the collision term of the transport equation. Each term in this expansion represents a different 
physical contribution to the transport equation. 



VI. EQUATIONS FOR A< 

Applying the translation rules (|28|l to Eqs. 122II23|I . neglecting the terms proportional to d^ due to the quasi-particle 
approximation (|27|l . and taking the difference and sum of the equations, we obtain 

'p^d^-\d,ns{x)d'A^/^^{x,p) = '-{n>{x,p)^<{x,p)-n<{x,p)^>{x,p)) 

- i{n^(X,p),A+(X,p) + A-(X,p)} 

- i{n+(x,p) + n-(x,p), A^(x,p)} , (30) 



m" - lisiX) A<{X,p) 



i (n^ [x, p) ( A+ (X, p) + A- (X, p)) + (n+ (X p) + n- (x, p)) a^ {x, p)) 
^{n>(x,p), A<(x,p)} - ^{n<(x,p), A>(x,p)} , 



(31) 



where we have used the identity H14|l applied to the Green's functions and self-energies. One recognizes Eq. (|3()|l as a 
transport equation while Eq. ()31|l is called a mass-shell equation. We note that in the limit of free fields, Eqs. (|3UI31f) 
become 



p^'d^A<{X,p) = Q 



(32) 



[/-m2]Ao^(X,p)-0 



(33) 



Due to Eq. H33() . Aq (X,p) is proportional to 5{p'^ — m?), and consequently free quasi-particles are always on mass-shell. 
We note that if the quasi-particle approximation H27|) were not used, the mass-shell equation would have the form 



rl 



-5^ -p^ +m' 



A§{X,p) = Q. 



(34) 



and the off-shell contribution to the Green's function An would be nonzero. 



The equations H3QI31|I can be rewritten in a more compact way: 

{p2_m2+n,(X) + Ren+(X,p), A^(X,p)} = t(n>iX,p)A<{X,p)-Il<{X,p)A>iX,p)) 

- [ll^{X,p),RcA+{X,p)}, (35) 

p'^~m^ + ns{X)+RcIl+{X,p)\A^{X,p) = -U^ {X , p)ReA+ {X , p) 

- l{u>{X,p), A<{X,p)] + l{u<{X,p), A>{X,p)] , (36) 



where 

ReA±(X,p)^l(A+(X,p) + A-(X,p)) , ImA±(X,p) ^ ±1(a+(X,p) - A-(X,p)) . (37) 

We note that (iA^{X,p)^ — —iA^{X,p) which is a direct consequence of the definitions (|11I12|I . The equation 
analogous to Eq. (|35|l was earher derived in [la . Il4l |25| . 

The gradient terms in the right-hand-sides of Eqs. H35I36|I are usuaUy neglected. The justification for dropping 
these terms is considered below. The small parameter that characterizes the gradient expansion will be denoted e. 
In Sec. IXIII wc will calculate the self-energies using a perturbative expansion in the coupling constant g. For the 
purposes of the discussion below, we will use a new parameter to characterize the coupling constant expansion. This 
new parameter is defined by \ = g for (p'^ and A = ^ for (p^ theory. The parameter A is introduced to simplify the 
notation: the leading order contributions to Ils{X) + Ren+(X,p) and Il^{X,p) are of order A^ and A'' respectively, 
in both (j)^ and (p^ theory ^. 

Now, we compute the order of each term in Eq. H35|l . We can split the Poisson-like bracket on the left hand side 
into two pieces. The first piece does not contain the self-energy and is of order e. The second piece is of order A^e. 
The first term on the right hand side (the term without the Poisson-like bracket) is of order A*, and the term with 
the Poisson-bracket (the gradient term) is of order A'^e. Thus, we find that the gradient term on the right hand side 
of Eq. H35|l is of higher order than the remaining terms and can be dropped. Similarly, the gradient term on the right 
hand side of Eq. (|36|1 is of order A^e and it is of higher order than the remaining terms. 

We note that the argument presented above breaks down close to equilibrium. At equilibrium the non-gradient term 
on the right hand side of Eq. (|35|l is identically zero, and thus cannot be considered to be bigger than the gradient term. 
In the equilibrium case, it has been shown that the gradient term does not need to be dropped but can be combined 
with the interaction term on the left hand side to produce the usual Vlasov term |2!tI |. The very-close-to-equilibrium 
situation is not fully understood. 

VII. EQUATIONS FOR A± 

We also write down the transport and mass-shell equations satisfied by the retarded and advanced Green's functions. 
Starting with Eqs. (|24I25|I . one finds 

|p2 _ ^2 _^ ^^(A) + n±(X,p), A±(X,p)} = , (38) 

p^ -m^ + Us (X) + n± {X, p)] A± {X, p)^l. (39) 

We observe that the gradient terms drop out entirely in Eq. (|39|l . The leading order corrections to this equation are 
second order in the gradient expansion. Equation H39|) can be immediately solved to give 

^* (^' P) = .^TT !y^^U±(Y ^ ■ (40) 

p^ - m^ + Us[Xj + U^[X,p) 

We note that any function f{K) satisfies the equation {K, f{K)} = 0, for K an arbitrary function of X and p. Thus, 
the expression igOI) solves Eq. ^^ as well as Eq. (p!?^ . 

The real and imaginary parts of A , which we will need later, are 

ReAHx,p) = P^-m^ + nsiX)+ReIi-iX,p) 

{p2 _ ^2 _^ n5(x) + Ren+(x,p)) + (imn+(x,p)) 

ImA±(X,p) = Tln.n^{X.p) ^^2) 

(p2 - to2 + ns{x) + Ren+(A:,p)) -i- (imn+(A:,p)) 

Using the relations (|11I12|I . we can obtain the Wigner transforms of the retarded and advanced functions for free 
fields 

A J {X, p) = / . (43) 

p'^ — m^ ± ipoO 

Comparing the expressions H4U|) and (|43|l . we find that in the limit g ^ the self-energies must satisfy 

ln.U-iX,p)^-lmU-iX,p)^{i: [^ ^^ > , (44) 



VIII. SPECTRAL FUNCTION 

In this section, we introduce another function which will be useful later on. The spectral function A is defined as 

A{x,y) =' mx),<j>{y)]) = iA>ix,y) - tA<{x,y) = iA+{x,y) - iA-{x,y) = T2ImA±(x,2/) , (45) 

where [(f>{x) , (j){y)] denotes the field commutator. For real fields, which are studied here, A{X,p) = —A{X,—p). 
From the transport and mass-shell equations (|35II36|I . one immediately finds equations for A{X,p) of the form 

{p2 _jri^ + Us{X) + ReU+{X,p), A(X,p)} = 2{lmU+{X,p), ReA+(X,p)} , (46) 

p^-m^ + lis {X) + Ren+ {X, p)] A{X, p) = 2 Imn+ {X, p) ReA+ {X, p) . (47) 

Substituting the formula 141() into the algebraic equation (|47|l . we find the solution 

A{X,p) = ^ '^' „ ^. 48 

(p2 _ ^2 + n^^x) + ReU+iX,p)y + {Imn+ {X,p)y 

It is easy to show that the function (|48|l solves Eq. H46|) as well. In fact, the spectral function (|48|l can be found 
directly from Eq 142|l due to the last equality in Eq. H45(l . 

The free spectral function can be obtained from the free retarded and advanced functions H43(l by means of the 
relation (|45|) . It can be also found from Eq. (|48|l using the condition (|44|l . In either case one needs to make use of the 
well-known identity 

P- T iT^S{x) . (49) 



x ± J0+ X 

We obtain 

Aa{X,p) = 2n6{p^ - m2)(e(po) - e(-po)) • (50) 

When Imll''' = the spectral function describes infinitely narrow quasi-particles with energies uniquely determined 
by their momenta. When Imn+ ^ the quasi-particles are of finite life time, and the spectral function gives the 
energy distribution of a quasi-particle with a given momentum. 

IX. DISTRIBUTION FUNCTION 

The distribution function f{X,p) is defined through the equation 

e{po) A{X,p) f{X,p) "^^ e{po)tA<{X,p) , (51) 

where A{X,p) is the spectral function H48|l . This definition is motivated by Eq. (|17|l which leads us to identify A> as 
the quantum analogue of the classical distribution function. We note that f(X,p) depends not on the three- vector p 
but on the four-vector p. Using Eq. H45|) and the identity (|15|l in the form A'^{X,p) — A^{X, — p), we have 

iA>{X,p) - e{pa) A{X,p) {fiX,p) + 1) - e(-po) AiX,p) fix, -p) , (52) 

iA<{x,p) = e(po) A{x,p) fix,p) - e(-po) A{x,p) {!{x,-p) + 1) . (53) 

In the case of free fields, Eqs. (|52I53|I simplify to 

iA>{X,p) = 2^5{p^~m'')(Q{po) (/o(X,p) + l) + e(-po) /o(^,-p)) (54) 

= ^5{Ep-po){UX,p) + 1) + ^5{Ep+p^)UX,-p) , 

iA<{X,p) = 27TS{p'-m')(e{po)fo{X,p) + e{-po){MX,-p) + l)) (55) 

= ^SiEp - po)foiX,p) + ^S{Ep + po){fo{X,-p) + 1) , 
rtip tLp 
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where the explicit form of the free spectral function (|50() was used. 

For future use, we write down the results for the propagators Ag, Ag and Aq^™ in terms of the free distribution 
function /q. Using the free retarded and advanced functions 1)43(1 and the identity 1(49(1 . we obtain 

zAg(X,p) = tA+{X,p) + tA<{X,p) 

= T--?^ + 2^ ^(P^ - ^^) (Q(po) /o(^' p) + Q(-Po) /o(^' -P)) (56) 

zAg(X,p) = -*Ao-(X,p) + *A<(X,p) 

= 2 , ^2', jo+ + 2^ ^(P' - "^') (©(Po) /o(^, P) + e(-po) fo{X, -p)) , (57) 

zA^^-(X,p) = *A>(X,p) + zA<(X,p) 

= 2^ V _„,2)|^2e(po) /o(X,p) + 2e(-po) foiX,-p) + 1) . (58) 

We note that in the case of equilibrium (homogeneous) systems, the familiar definition of the distribution function 

is 

r»1^'i(at(p)a(p)), (59) 

where V denotes the system volume and a(p) and a^(p) are annihilation and creation operators, respectively. The 
definition l|^ should reduce to l|5^ when an equilibrium (homogeneous) system is considered. We discuss this point 
below for the case of a noninteracting system for which the field that solves the equation of motion ||2J) can be written 
as 

where k = (Ek , k) . Substituting the solution l(60() into the definition (jB)) one finds 

tA>{X,p) ^ I /""^ -^ '^''^ {2iTf (61) 

"^ '"^^ J ^i2n)^2Ek ^{2^f2E, ^ ' ^ ' 



^(^)^/=f:(^"^^^"(^) + ^^^'^"^^^))' («o) 



-»(fe+,)X^(4) (^p _ fc_^^ (a(k)a(q)) + e-»('=-?)^5(4) {p - ^) (a(k)at (q)) 

+e'(fe-.)X5(4) (^p ^ fc + 9^ (at(k)a(q)) + e^('=+'^)^5(4) [p + ^—1^ (at(k)at(q)) 

Since we restrict ourselves to the consideration of homogeneous systems, we require that the function A^{X,p) is 
independent of X. One ansatz that satisfies this constraint is given by 

(a(k)a(q)> = , (at(k)at(q)) = , (62) 

and 

(a(k)at(q)) = ^_^(a(k)at(k)) , (ank)a(q)) = ^^(at(k)a(k)) . (63) 

Using these expressions, the integrals in Eq. ((61(1 are trivially performed and we obtain 

zA>(p) = :^5{E,-p,) i(a(p)at(p)) +^5{E,+po) i(at(-p)a(-p)) . (64) 

Using the commutation relation [a(p), a'l'(k)] = (5^^-'(p — k), which gives 

(a(k)at(k)) = (at(k)a(k))+F, 

one finds that Eq. ((64(1 combined with the definition 1(59(1 reproduces ((5511 . Thus we have shown that the definitions 
((51(1 and ((59(1 are consistent with each other for homogeneous systems. One notes that for nonhomogeneous systems 
the conditions I(62I63() are, in principle, not fulfilled. As a result, the functions A^{X,p) are not only X— dependent 
but also have support for off-shell momenta p in agreement with Eq. I(34() . Therefore, the expressions ((54155(1 can be 
treated as a representation of Aq(X,p) only for weakly nonhomogeneous systems. 
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X. MASS-SHELL CONSTRAINT 



In this section we discuss two issues that are related to the mass sheh constraint. We first show that in the 
homogeneous hmit (or neglecting gradient terms), the mass shell condition is satisfied if the distribution function 
is defined in terms of the spectral function as in Eq. H51|) . and if this distribution function satisfies the transport 
equation. Secondly, we show that using this definition of the distribution function, the mass-shell constraint reduces 
to the familiar definition of the mass-shell condition on the four-momentum in the limit of zero width quasi-particles. 

The transport and mass-shell equations for a homogeneous system, which are obtained from Eqs. 1)35(1 and ((36|l by 
dropping gradient terms, read 



n>{x,p)A<{x,p) - n<{x,p)A>{x,p) 

p^-m^ + Us {X) + Ren+ {X, p)]a^ {X, p) 







-n<(X,p)RcA+(X,p) 



(65) 
(66) 



These equations can be rewritten using Eqs. (|52I53|) which express A'^{X,p) through the distribution function. We 
work below with the mass-shell equation for A<(X, p) and restrict to po > 0. It is straightforward to repeat all steps 
with po < 0, and for A< with po > and po < 0- We find 



A{x,p) n>(x,p)/(x,p)-n<(x,p)(/(x,p) + i) = o, 

p^-m^ + I].s{X) + Y{.cI].+ {X,p)\A{X,p)f{X,p) = -tIl<{X,p)ReA+{X,p) . 
Using the spectral function equation H47|) . the mass-shell equation (|68|) has the form 



ReA+(a;,p) 



n>ix,p) fix,p) - u<ix,p) {fix,p) + 1 



= 



(67) 
(68) 

(69) 



One sees that if / solves Eq. H67|l . it automatically satisfies Eq. ((69|l . Equivalently, the transport and mass-shell 
equations (|35|l and (|36|) can be replaced with the two equations H47|l and (|67|l in the homogeneous limit. 

We note that the equations H67I I68|l . which were obtained by neglecting gradient terms, are exact for equilibrium 
(homogeneous) systems. These equations can be written as 



n>(p)r'>(p)-n<(p)(r^(p) 



(70) 



which is the well-known Kubo-Martin-Schwinger (KMS) condition, see e.g. jlSl], that is satisfied by the self-energy in 
equilibrium. We also note that in going from equations l|5S|l and l|^ to equations l|17|l and (|^ . we have replaced four 
equations with two. This reduction in the number of independent Green's functions is expected since in equilibrium 
the two functions A^{X,p) are not independent but related through a KMS condition of the form ((70|l . 

We have shown above that when the distribution function is defined by Eq. H51|l . the transport and mass-shell 
equations H35I36|) can be replaced by the equations (|47|) and 16 7|) . In the limit Imn+ -^ the solution to Eq. H47|l . 



(71) 



which is given by Eq. (|48|l . becomes 

AiX,p) = 2^(e(po) - e(-po)) Sip^ -m^+ HsiX) + RcU+{X,p)) . 
The argument of the delta function gives the usual form of the mass-shell constraint of the four-momentum 

p2 - m^ + Ils{X) + RcIl+{X,p) = , 



(72) 



which tell us that only three out of four momentum components are independent. For Imn"*" small but finite, this 
equation determines the position of the maximum of the spectral function. 

We note that when the system of interest is significantly inhomogeneous, cf. Eq. H34|) . or when an interaction 
generates a non-negligible value of Imn"*", the function A^{X,p), which solves Eq. H36|l . has support from momenta 
not satisfying the relation (|72|) . In this case, one has finite-width quasi-particles and the equation H72|l only gives, 
according to Eq. H48|l . the most probable energy of a quasi-particle. Thus, the statement that the four- momentum is 
on the mass-shell (|72|l and the statement that the function A<(X,p) satisfies the mass-shell constraint (|36f) are, in 
general, not equivalent. 
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XI. TRANSPORT EQUATION 



The distribution function / satisfies the transport equation which can be obtained from Eq. H35|) for A> or A<. 
Using Eqs. H46|) and (|51|l . one finds 

e{po)A{X,p){p^^m^ + nsiX) + Ren+ {X,p), f{X,p)} 

= iQ{po) [a{x,p) (n>(x,p) fix,p) - ii<ix,p) {fix,p) + 1)) 

+ f{X,p) {n>(X,p), RcA+(X,p)} 

- (/(X,p) + 1) [n<{X,p), ReA+iX,p)}] , (73) 

where we have used the identity 

{A, BC} = {A, B}C+{A, C}B . 
Neglecting the gradient terms in the r.h.s. of Eq. 173(1 . we obtain the equation 

e(po){p' - m.^ + Usix) + Ren+(x,p), f{x,p)} = ie(po)(n>(Xp) fix,p) - n<{x,p) {!{x,p) + 1)) . (74) 

The transport equation H74|l greatly simplifies for zero- width quasi-particles. Using the mass-shell condition 172|l . 
and taking the positive energy solution, one can evaluate the Poisson-like bracket to obtain 

e(po) [p" - m^ + ns{X) + Rcn+(X,p), /(X, p)} ^ 2£;p(^ + v V)/(X, p) + 2VV{X,p) Vp/(X, p) , (75) 

where 

V{X,p) = I{s{X) + RcIi+{X,p), 



and the velocity v = dEp/dp. Note that Ep is not -^/p^ + m? but the positive solution of Eq. H72|l . Substituting 
Eq. 175|l into Eq. ((74|l . the transport equation takes the form 

i^p(^ + v-v)/(x,p) + vF(^,p)Vp/(x,p) = e(po)^[n>(x,p)/(x,p)-n<(x,p)(/(x,p) + i)] . (76) 

In the left hand side of Eq. H7()|) . one recognizes the standard drift and Vlasov terms [33 ■ The calculation of the 
self-energies Ren+ and XI^ in the Vlasov term is straightforward and the physical interpretation is well known: they 
contribute to mean field effects and result in the generation of an effective mass. In this paper we will not discuss 
these terms (they have been studied in detail in 18]). From the right hand side of Eq. H76|l . one defines the collision 
term C{X,p) as 

e(po)c(x,p)1^fe(po)^[n>(x,p)/(x,p)-n<(x,p)(/(x,p) + i)] . (77) 

The structure of this collision term is the subject of the next five sections. 

XII. COLLISION TERM STRUCTURE 

In order to obtain physical results from the transport equation derived in the previous sections, we must specify the 
self-energies that appear in the collision term. We will assume that the interaction is sufficiently weak that the self- 
energies can be expanded in powers of the coupling constant, and thus expressed through the Green's functions. We 
will develop a diagrammatic expansion of the self-energies and identify the different physical processes that contribute 
to the collision term (|77|l . The collision terms we obtain are local in space-time, and consequently valid in the 
Markovian limit where all memory effects are neglected. 

We first rewrite Eq. lf77|) by multiplying by the spectral function A{X,p) and using Eq. ((^ . Thus, we obtain the 
equation 

e(po)A(x,p)c(x,p) = e(po)^(n<(x,p)A>(x,p)-n>(x,p)A<(x,p)) . (78) 
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Evolution along the contour illustrated in Fig. ^ is formally very similar to evolution along the real time axis, 
and consequently it is possible to define perturbation theory on the contour. Self-energies with real time arguments, 
in particular n<, have to be extracted from the contour self-energy. Since these calculations involve summations 
over contributions from both branches of the contour, they are much more difficult than their vacuum counterparts. 
There are several different methods that can be used which are based on different representations of real time field 
theory |23| and each has its advantages and proponents. At the lower loop levels, the relevant self-energies have been 
calculated in several different ways both in equilibrium 39] and non-equilibrium |l8||. (A general structure of the 
equilibrium self-energy expressed through the scattering amplitude has been studied in |40|.) At three or more loops, 
the calculation is extremely cumbersome and requires the use of special methods. 

Our calculation of the self-energies, which enter the collision term of transport equation, uses the Keldysh repre- 
sentation formulated in terms of retarded, advanced and symmetric Green functions. To understand why the Keldysh 
representation is better suited to perform such a calculation, we need to look at the Wigner transformed propagators 
for non-interacting fields in the Keldysh and 1/2 representation. In the 1/2 basis the propagators are given by equa- 
tions (|5^ l55l l5?n [KTjl and in the Keldysh basis by (031 EHJ- The functions iAg and iAq do not contain off-mass-shell 
contributions that have the form of propagators. They are non-zero only on the mass-shell. In contrast, lA^ and iAq 
contain one piece that corresponds to time-ordered and anti-time-ordered propagation, respectively, and one piece 
that is non-zero only on the mass shell. In the Keldysh basis, however, none of the Green's functions mixes the 
propagating pieces and terms that are non-zero only on the mass-shell: Aj and A(^ contain only the contributions 
corresponding to retarded and advanced propagation, respectively, and Aq^™ is non-zero only on-mass-shell. Remem- 
bering that on-mass-shell Green functions represent real particles and propagating Green functions virtual particles, 
it is not unexpected that this clean separation of real and virtual particles makes the physical interpretation of the 
collision term much more straightforward. 

We switch to the Keldysh basis using the relations H14|) and analogous relations for the self-energies. The collision 
term (|78|) . which will be computed in the next sections up to four- loop level, then equals 

eipo)Aip) c{p) - ie(po) [ff^'"b)(A+(p) - A-(p)) - (n+(p) - n-(p)) A^y-(p)] . (79) 

It will be shown that in the Keldysh representation, all contributions to the collision term have the form of 'cuts' 
[27j . We use the word 'cut' in the way that it is usually used when discussing the calculation of the imaginary part 
of a diagram: a 'cut line' passes through a diagram from top to bottom dividing it into two parts, so that each part 
contains at least one external leg. A propagator that is crossed by the cut line becomes a 'cut propagator' which 
means that it has been put on the mass shell. 

We introduce the notion 'central cut' to refer to the cut for which all loops contain at least one cut propagator. 
For the one-loop diagram in cf)^ theory, and the two-loop diagram in 0^ theory, the only possible cut is the central 
cut. At two and more loops in (j)^ theory, and three and more loops in 0^ theory, all diagrams contain contributions 
from non-central cuts. In addition, there are some diagrams for which no central cut exists. In each case however, 
the contributions from non-central cuts can be included as renormalized lower loop terms. An example of this is 
presented in Sec. lXIVBl The main thrust of our paper involves the analysis of the central cuts of three- and four- loop 
diagrams. 

The calculation of a self-energy graph involves two non-trivial steps: the calculation of the integrand and the 
computation of the integral itself. In this paper we are primarily interested in the calculation of integrands which 
determine the structure of the collision term. At three or more loops, calculating by hand in any representation is 
extremely tedious. We use a MATHEMATICA program developed by one of us, which is described in detail in [26], to 
calculate the integrand. To use the program one assigns momenta and indices to every line and vertex in the diagram. 
Each index can take the values 1 or 2, corresponding to the top and bottom branches of the contour. The indices 
that correspond to external legs take fixed values, and a specific combination corresponds to a specific self-energy in 
the 1/2 representation, in analogy with Eq. ^. Internal indices are summed over their two possible values 1 and 2. 
Using Eqs. H8I13I14|I . one switches to the Keldysh representation. A huge number of terms is produced even for fairly 
simple diagrams. However, there are many cancellations between these terms. The program identifies the surviving 
contributions and provides them as output. 

In Sees. IXIlTl and lXIVl we rederive the well known one- and two-loop results to illustrate our technique. In Sec. IX VI 
we use our method to obtain results for three-loop diagrams which contribute to 2 <-^- 3 processes in (jy^ theory, and in 
Sec. lXVll we study four-loop diagrams in t/)^ theory which correspond to 2 <^ 4 and 3 <-» 3 processes. 
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XIII. ONE-LOOP CONTRIBUTIONS 

We start by considering the lowest order contributions. In the 0^ model, the only one-loop diagram is the tadpole. 
It is well known that this diagram only contributes to mean-field dynamics, see e.g. |18| . and therefore it will not be 
discussed here. For (j)^ theory, the one-loop contribution to the self-energy, which is shown in Fig. [3 has a non-trivial 
imaginary part. The related scattering amplitude corresponds to 1 <-!■ 2 processes which are kinematically forbidden 
for on-mass-shell particles, as an on-shell particle of mass m cannot decay into two on-shell particles with the same 
masses. However, it is important to remember that the process is allowed for virtual particles and has been repeatedly 
discussed, in particular, in the context of electron-positron-photon interactions. We will discuss this contribution in 
detail to illustrate the method that we will use in the analysis of more complicated multi-loop diagrams. 

At this point, we introduce some notational simplifications. We suppress the X— dependence of all functions. All 
Green's functions used in our perturbative calculations correspond to free fields, and thus we suppress the index '0' on 
these functions. Similarly, we suppress from now on the index '0' which has been used to denote the free distribution 
function and the free spectral function. Finally, we define on-shell four-momenta 

P=(^P,P), (80) 



with Ep = \/p2 -t- w?. We note that \/p-^ -f vn? should be denoted £?, but as above the index '0' is suppressed. 
The Keldysh self-energies corresponding to the diagram from Fig. |2l are 

n+b) = *Y / ^[A-(9)A^^'"(9 +P) + A^>'-(g)A+(g + p)] , (81) 

n-(p)-*^y'^[A^'""(g)A-(g+p) + A+(g)A^y-(g+p)], 

ff'""(p) =«yy'^[A^''"X9)A^^'"(g + p) + A+(q)A-(g + p) + A-(p)A+(g + p)] . 

We substitute these self-energies into the right hand side of the collision term H79I) and express A''y™(p) through the 
distribution function according to Eq. (|58|) . We also use Eqs. H45II5U|) which give 

z(A+(p) - A-(p)) = sgn(po)2^,5(p2 _ m") , (82) 

and we observe that 

'^''^ A+(<z)A+(g+p)= f-^A-{q)A-{q+p)=Q, (83) 



(27r)4 '^' ^^ ' ^^ ; (2^^4 
due to the positions of poles of A^. Thus, we obtain 



a f da 

e(po) A(p)CiL-03(p) = e(po) Y^ / 7^-j4 X! •^3(n-p,P;n.5,q;-rir,-r) (84) 

X If tIq sgn(ao) 27r5(a^ — m^) 8(nQao) , 

a— {p, q. r] 

where r = q + p and a is a generic four-momentum variable which equals p, q or r, and Ua = ±1- The statistical 
factor J^N is defined as 

J^Af(ni,pi;n2,P2;- • -nwjPJv) = fl{ni,pi)n{n2,P2) ■ ■ -^{nN^PN) (85) 

- il{-ni, -pi) V,{-n2, -P2) • • • ^{-UN, -Pn) 

with 

n{n, p) = (1 + n) (1 + /(p)) + (1 - n) /(-p) . (86) 

We note that the asymmetry in the signs of {ug, q) and {ur, r) in the factor Tz{np, p ; rig, q ; — n^, — r) occurs because 
the momenta q and r are defined to flow in opposite directions in Fig. [3 




q+p 
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FIG. 2: One-loop contribution to the 4> self-energy. 



To further simplify our result, we write the delta functions as 

sgn(ao) 5{a'^ - m^) = — — ^ rv e(rvao) ^("o - riaEa) , 



2Ea 



(87) 



n„=±l 



where Ea = ^/m? -f a^, and get 

e(po)A(p)Cii_^3(p) = e(po) 



5^ 
16 



(2^ 



y^ J^3(np,p;n5,q;-nr,-r) ]^ — 5 (ao -«„£'„) (88) 



i?Q 



{np,nq,nr}=±l a={p,q,r} 

We rewrite the expression H88|l by introducing the delta function S^'^' {j) + q — r) and an integral over r, which allows 
us to treat r as an independent variable. Performing the integrals over po and tq, one obtains 



e(po) A{p) C,,^,sip) = Qipo) i^jA3_jAl_ (2^)3^(3) (p + q _ r) [] ^ 



(89) 



Q — {p, g, r} 



X X! 5(po - 'T-p-E'p) 27r5(npi?p +n,£:q - n^i?,-) .F3(np,p;ng,q;-nr,-r) 

{np,nq,nT.} — ±1 

We discuss below the physical interpretation of the positive energy piece of the collision term (|89f) which will be 
denoted with the index (-I-). We drop the theta function in Eq. 189|) and take only the Up = 1 term in the sum. We 
also use the explicit form of the free spectral function (|50ll to make the replacement 

eipo) A(p) ^ ^ S{po - Ep) . (90) 

In principle, we are left with four terms that come from the sum over n^ and rir in Eq. H89I) . One of these terms 
corresponds to a delta function of the form S{Ep + Eq + E^) which has no support. The remaining three terms 
correspond to the processes symbolically denoted as: p^^q + r^p + q-i^r and p + r '^ q. The last two terms can be 
combined by making the variable transformation q <-> r in the second term. 



We consider first the contribution to the collision term H89|) obtained from choosing 



-1 and Ur — 1. Making 



the change of variables q ^ — q, we find 

.2 , ^3q 



C. 



\ T J.JP ^-^ <? + ''I = — 



(Pr 



{2Tr)H^^\p-q-f)T^{\,V--l,-ci;-l,^v) 



{2T:fEq J (27r)3£;, 
where the statistical factor equals 

.^3(1, p ; -1, -q ; -1, -r) - 2^ [(l + /(p))/(q)/(r) - /(p) (l + /(q)) (l + /(r)) 

This process is represented in Fig. |31 



(91) 

(92) 




FIG. 3: p ^^ q + r scattering process. 
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FIG. 4: p + q ^^ r scattering process. 

Next, we look at the contributions to Eq. H89(l obtained from choosing Uq = rir = 1 and Uq = Ur = —1. After 
changing variables and combining the two terms, we get 



with 



CJt^3b + .^.]^f;/^/^(2.r^(^)(p + .-.1^3(M;i,^ 



.F3((l, p ; 1, q ; -1, -r) - 23 (1 + /(p)) (l + /(q))/(r) - /(p)/(q) (l + /(r)) 



This contribution is illustrated in Fig. 0] 
The complete one-loop collision term equals 



C. 



(+) 

IL-th^ 



(P) 



d^q 



dh 



i2Tr)3Eq J {2n)^Er 



\M\'s 



(93) 



(94) 



(95) 



(2^)45W(p - q - r) - [(1 + /(p))/(q)/(r) - /(p) (l + /(q)) (l + /(r)) 



2! 



+ (2^)V4)(p + g-r~) [(l + /(p))(l + /(q))/(r)-/(p)/(q)(l + /(r)) 



where the scattering matrix element equals |Al|ii_^.'i = 5^/4. The two terms proportional to /(p) correspond to 'loss' 
contributions and the terms proportional to (1 + /(p)) correspond to 'gain' contributions. In the first term, which 
represents the process p ^^ Q + f, there is the factor 1/2! because of the integration over momenta of two identical 
particles in either final or initial scattering state. 

We remind the reader that the contribution (|95|l to the collision term is identically zero because it is impossible 
to satisfy the delta function constraint. Physically, an on-shell particle of mass m cannot decay into two on-shell 
particles with the same masses. We have presented this calculation to illustrate the method that we will use in the 
analysis of more complicated multi-loop diagrams. 

XIV. TWO-LOOP CONTRIBUTIONS 

First we note that at any number of loops greater than one, there are two different kinds of diagrams that contribute 
to the self-energy. The first class contains diagrams that are built from diagrams of lower loop order by adding tadpole 
type insertions. An example is shown in Fig. |S| The contribution to the collision term from a diagram of this type 
has the same structure as the corresponding lower loop diagram with the addition of effective masses coming from the 
tadpole insertions. This type of diagram will not be discussed here. The second class of diagrams does not contain 
tadpoles. This is the type of diagram that we study in this paper. For the first few diagrams we give some of the 
details of the calculation but for the most part we focus on the final results. 




FIG. 5: Diagram with a tadpole insertion. 
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At the two-loop level, there are two contributions to the self-energy in (f)^ theory and one contribution in (j)^ theory 
(excluding diagrams with tadpoles). 



A. (f)'^ theory 

The two-loop contribution to self-energy is shown in Fig.|Sl As in the case of the one- loop contribution in cjr' theory, 
there is only one cut which is the central one. Calculating the retarded, advanced, and symmetric components of the 
self-energy, we obtain 



n+(p) 



9 



d'^q 



24 J (27r)4 J (27r) 



dH 



(96) 



n-W-fr 



dS 



24 J (27r)4 J (27r 



dH 



(97) 



A-(g)(A"y'"(/)A"y'"(s) + A'{l)A+{s)) 
A^y'"(g)(A-(/)A^y'"(s) + A^y™(OA+(s)) 
A''y™(s)(A''y"(g)A+(0 + A'^y'"(OA+(g)) 
A~(s)(A"y"^(OA"y"^(q) + A+(/)A+(q)) 

A+(/)(A^(s)A^>'"^(g) + A^y'"(s)A+(g)) 

+ A-(0(A-((?)A^y'"(s) + A^>""((?)A+(s)) 

+ A''y™(0(A''^™(g)A"y'"(s) + A'{s)A+{q) + A-{q)A+{s)) 

where we have used s = p + q + l. We substitute these formulas into the collision term (|79() and express the symmetric 
Green's function through the distribution function according to Eq. H58|l. As in the case of the one- loop diagram, we 
use Eqs. 182I83II to find 



u^y-^ip) = ^ 



d^g 



24 J (27r)4 J (27r 



dH 



e{po) Aip) C2L-ripo) = e(po) 



5V 
192 



d^q 

(2^ 



dH 
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(99) 



y^ T4(np,p;ng,q;ni,\;-ns,s) JJ — S{ao - UaEa 



,ni ,ns } — il 



a={p,q,ls} 



Ea 




FIG. 6: Two-loop contribution to the (f'^ self-energy. 

We rewrite the result (|99() in a more symmetric way by introducing the delta function 5^^^ {p + q -\- 1 — s) and an 
integral over s, which allows us to treat s as an independent variable. As in the previous section, we use the free 
spectral function H9U|I and consider only positive energy contributions to the collision term. Performing the integrals 
over go, 'oi and sq, we have 



c'Ci^-Ap) 



5%- 1 



192 (2 



L_ J d^q J dH J dhS^'Hp + n-i-s)^ 



qElEg 



(100) 
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X 22 5{Ep + nqEq + niEi -UsEs) T4il,p;ng,q;ni,l;-ns,-s) 

{np,ni.ns}=±l 

In principle the sum produces eight terms, but five of these eight terms are zero because the corresponding energy 
conserving delta function has no support. There are three choices of {nq,ni,ns} that give non-zero contributions 
to the collision integral. These three terms can be combined by making an appropriate variable transformation on 
each term. We list below: the three choices for the Ui variables with i = {q, Z, s} that give non-zero contributions, 
the variable transformations for each of these three terms that allow us to combine them, and the corresponding 
transformations to on-shell four- momenta (using the notation H8U|) '): 



[1] Hq = -ni = Us = I; 
[2] Uq = -ni = -ris = -1; 
[3] Uq = ni = lis = -1; q 



1-^-1; 




Q- 


^ q, I ^ 


-I, s- 


■* s; 


q^-1, l^q; 




q- 


.-1,1^ 


-^ q, s ~ 


-> s; 


-s, 1^-1, s^ 


q; 


q^ 


-S, I -^ 


-I, s- 


-> — 



(101) 



The statistical factor and the two delta functions have the same form for each term. Combining the deltas into 
one four-dimensional delta function that expresses energy-momentum conservation and using the explicit form of the 
statistical factor defined by Eqs. (|85I86|I . we finally find 



a 



(+) 

2L- 



,^[p + q^l + s] = - 



d\ 



dH 



dh 



(27r)4,5(4)(g+p-f-S)|A^| 



2 
2L- 



(102) 



1 

2! J {2Tr)^Eq J {2Trf El J {2tt)^E, 

1 + /(P)) (1 + /(q))/(l) /(s) - /(P) /(q)(l + /(I)) (1 + /(s)) 

i2^_a,4 — 5^/16. Eq. (|102() gives the standard collision term corresponding to the binary process symbolically 
denoted asp + q^^l + s and illustrated in Fig. [7| The first (second) contribution is the usual gain (loss) term. 



where \M\? 




FIG. 7; p + q ^^ I + s scattering process. 



B. (j>^ theory 
In (p^ theory, there are two graphs of different topologies that contribute at two-loop order as shown in Fig. |S1 





FIG. 8: Contributions to the self-energy in (jf theory. 

We begin with a detailed calculation of the contribution to the collision term from the first diagram in Fig. ^2 
The results for the self-energies are as follows 



n+(p)-^ 



d'^q 



d'^k 



(27r)4 J (27r) 



A-(q)A~(r)A'^y"^(r)(A-(/)A"y'"(fc) -I- A''y'"(OA+(fc)) 
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(103) 
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(104) 



nsy-(p) ^ ^ 



d^g 



d^fc 



(27r)4 y (27r) 



(A-(r))^A"y"X9)(A"(0A'^"'(fc) + A'^>""(/)A+(fc)) 

+ A+(r)A^(r)A''y'"(r)(A-(/)A"y™(fc) + A''>''"(OA+(fc)) 

+ A-(r)A+(p)A+(r)(A'^y"Xfc)A'^'"(0 + A-(OA+(fc) + A-(fc)A+(/)) 

+ A''y'"(r)A+(p)A+(r)(A-(/(;)A'^y'"(0 + A">""(fc)A+(0) 

A-{r)A'y"'{r)A'y'^{q){A-{l)A'y'^{k) + A"y'"(/)A+(fc)) (105) 

+ A+(r)A"y™(r)A''y'"(q)(A~(fc)A">'™(/) + A'^y'"(fc)A+(0) 

+ (A-(r))'A+(g)(A-(OA''y''^(fc) + A^y'°(/)A+(fc)) 

+ (A+(r))'A-((7)(A-(fc)A^>'"^(/) + A^>""(fc)A+(/)) 

+ A-(r)A+(r)A''>""((?)(A">'"^(fc)A">'"XO + A-(/)A+(fc) + A-(fc)A+(0) 

We substitute these expressions into the colhsion term (|79|l and divide the resuhs into contributions of three different 
types, which correspond to contributions to the three cuts shown in Fig. O 




(a) 



(b) 



(c) 



FIG. 9: Cuts of the first diagram in Fig.|Hl 

These cuts are referred to as the 'central', 'left' and 'right' cuts, respectively. Terms that contain a factor A+(r) A^ (r) 
correspond to the central cut. The terms with A+(r) but without A~(r) are associated with the right cut, and terms 
that contain A^(r) but not A+(r) correspond to the left cut. The results read: 



e(po) A{p) c{p) 
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(108) 



^Q 



In Eqs. p07|) and (|108|l the factors 11+ (r) and n^(r) refer to one-loop self-energies of the form shown in Fig.|21 Thus, 
we find that the non-central cuts give contributions that can be understood as corrections to the one-loop result as 
shown in Fig. ^| 
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FIG. 10: Non-central cut for the first diagram in Fig. |H1 

Non-central cuts can always be treated in this manner i.e. every non-central cut provides a higher order correction 
to the central cut of some other diagram. As explained in the Introduction, we are only interested in the dominant 
contributions to any given physical process, and consequently we consider only central cuts from now on. We calculate 
contributions to the collision term at leading order in g for each physical process from the relevant sum of centrally 
cut self-energy diagrams at the same order of 5. The collision term contains some phase space integrals, the statistical 
factor JF/v (Eq. (|85(l ^. and the square of a matrix clement which is constructed from uncut propagators A . At higher 
orders in g, there are many self-energy diagrams that have to be included, some of which have more than one central 
cut. It is a non-trivial technical problem to combine all of these terms into a result which can be written as the square 
of a matrix element. This structure is obtained by an appropriate choice of momentum labels for the internal lines of 
the self-energy diagrams. We proceed as follows: 

1) the same set of momentum variables is always assigned to the cut propagators (equivalently, each central cut 
has the same !F factor); 

2) all possible permutations of momentum variables for the uncut lines are considered; 

3) for each diagram, all permutations are summed and the sum is normalized with the appropriate weight. 

Using this strategy, we calculate the nine diagrams shown in Figs. 1111 and 1121 where we have made the following 
definitions for momentum variables: 



r = p 



t — q + I , h — p + I , s = r + I 



(109) 



We note that the second topology in Fig. |S1 appears twice as often as the first topology (six times as compared with 
three) because of the fact that we take only one of two possible diagonal cuts for these diagrams. 




FIG. 11: Three permutations of the centrally cut first diagram in Fig. |H1 




21 




(d) 



(e) 



(f) 



FIG. 12: Six permutations of the centrally cut second diagram in Fig. |H| 



Defining the operation 

< • • • >2L- 



4 4 
dof „ . ^ g TT 



d*q 



dH 



E 

{np,ng,ni,ns}— ±1 



(27r)4y (27r)4 
J'4(np,p;n,,q;n;,l;-ns,-s) ]J — (5(ao - »^a-Ba) • • • 



a— {p,g,^,s} 

the contributions corresponding to the diagrams from Figs. lllM12t' can be written as 

CIt) = « A+(r)A-(r) »2l-03 , (Et) = « A+(/i)A-(/i) »2l_03 , dt) = « A+(t)A-(<) »2i_^ 
(Ell) = « A+(r)A+(i) »2L_03 , (d) = « A+(r)A-(/^) »2l_03 , (O) = « A+(/i)A+(t) »2l-^ 
df) = « A+(/i)A-(r) »2L-03 , d) = « A-(i)A-(r) »2l-03 , CUf) = « A-(t)A-(/^) »2l-^ 

After summing these contributions we obtain, 

e(po) Ab) C2L-03(p) = i « |A+(r) + A+(;i) + A-(i)|' »2L-03 . 



(110) 



(111) 



In order to express the result (|lll|l in a more symmetric way, we introduce the delta function S^'^' {p + q + I — s) 
and an integration over s, which allows us to treat s as an independent variable. As we have done previously, we use 
Eq. (|90|1 and consider only positive energy contributions to the collision term. Performing the integrals over pq, Iq 
and So J we obtain 



cSl,>(,» = "'-' 



l5/A/«'i/<i'.^™(c + P + i-.)j^ 



192 (2 



A+(p + q) + A+{p + l) + A+(-g - 



(112) 



{nq,?7;,ns}— ±1 



po = Eg, qo=ngEg, lQ=niEi 

5{Eq +npEp + niEi - UgEs) J^4(1,P ; n,, q; n/,1; -n^, -s) 



Except for the matrix element, the right-hand-side of this expression is the same as the right-hand-side of Eq. (|l()()(l . 
Eqs. H 1001 11 2(1 have similar structure since both correspond to contributions to the collision term from binary processes. 
As in the case of Eq. H100|) . the sum produces three non-zero contributions. The three choices of the n's and the 
variable transformations that allow us to combine terms are the same as in the 0* case and are given in Eq. (|101|l . 
The delta functions and statistical factors are the same as well. A straightforward but tedious calculation shows that 
the matrix element has the same form for each term. Thus, we finally obtain 
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The matrix element squared is given by 



A^ 



9 



1 1 1 

- + - + - 



where S, T, and U are the Mandelstam variables defined as 

S={q + p)\ T={l-pf, U=(q-lf 



(114) 



(115) 



We have used here capital letters to avoid confusion with momentum variables. This matrix clement is illustrated in 

Fig. 1121 
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FIG. 13; Scattering diagrams that contribute to the process p + q ^^ I + s \n (fP theory. 

XV. THREE-LOOP CONTRIBUTIONS 

A. (/)'* theory 

In (/)'* theory, the only three-loop diagram that does not involve tadpoles is shown in Fig. 114b .. This diagram does 
not have a central cut. The non-central cut presented in Fig. 114b . produces a correction to the two- loop sunset graph, 
as shown in Fig. [Tib , and will not be discussed. 




FIG. 14: Three-loop contribution to the (j)^ self-energy (a); and an effective two-loop contribution containing a corrected vertex 
function (b). 



B. 0^ theory 

There are diagrams of eight different topologies which contribute to the (f)^ self-energy at the three-loop level. They 
are shown in Fig. 1151 
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FIG. 15: Three-loop contributions to the (/>^ self-energy. 

We follow the same strategy as in the previous section and consider only the central cuts of the diagrams in Fig. ITKl 
The internal lines are labeled in all possible ways that produce the same statistical factor. In addition to the definition 
(I109|l . we define the following momentum variables: 



y — q + k , m = k + I , j — k + p , z = p + q + k , 
w = q + k + I , V = p + k + I , u = p + q + k + l. 
In analogy to Eq. H11U|) , we define the operation 
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Ec 



5{aQ - UaEa) 



{np,nq,ni,nk.n^}=±l 

Combining the results for the diagrams in Figs. [T5k - ll5h . we obtain 

eipo) A{p) C,L^^.ip) = l«|A+(t)A-(s)-f A+(y)A-(z) + A+(m)A-(t;) 

+ A+{m)A-{r) + A-{r)A-{z) + A-{r)A-{s) + A-{h)A-{v) 
+ A+iy)A-{h) + A-{h)A-{s) + A+{t)A'{j) + A-(z;)A-(j) 

+ A-{z)A-iJ) + A+im)A+iw)+A+it)A+iw) + A+iy)A+iw) »3l-03 



(118) 
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We can rewrite the result (|118() by introducing the delta function S^^^ {p + q + k + l — u) and an integral over u. As we 
have done previously, we use Eq. H9U|I and consider only positive energy contributions to the collision term. The sum 
over Uq^ ni, nk, and n„ produces 16 terms. Of these 16 terms, 10 are kinematically allowed. Four of these processes 
can be combined by relabeling momentum variables to produce the process symbolically denoted asp + g <-)■ l + k + u. 
The remaining six can be combined to produce the process of the form p + q + l'^k-'ru. 

The process p + q*r^ l + k + u 



We list the four possible choices of the n's and the variable transformations that allow us to combine terms and 
obtain the matrix element corresponding to the process p + q'^l + k + u: 



-ni 



-nk = n„ = 1; 1 
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[1] nq 

[2] Uq = ni = -Uk = -Uu = -1; 

[3] Uq ^ -ni =nk = -riu = -1; 
[4] Hq = ni = Uk ^ nu = -1; 



The delta functions and statistical factor have the same form for each of these four terms. A straightforward but very 
tedious calculation verifies that the same matrix element is produced in each case. The collision term representing 
the process p + q^^l + k + u equals 
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where the matrix element squared of the diagrams shown in Fig. 1161 is 
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A-{p + q) [A-(fc + l) + ^-{k + u) + A- {I + u)] 

+ A-{p- I) [A+{q -k) + A+{q -u) + A-{k + it)] 
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FIG. 16: Scattering diagrams that contribute to the process p + q^^l + k + u in the order they appear in Eq. 1121^ . 



The process q + p + l ^ k + u 

We list the six possible choices of the n's and the variable transformations that allow us to combine terms and 
obtain the matrix element corresponding to the process p + q + l'r^k + u: 



[1] Uq == ni = -Uk = riu = 1; k -> -k; 

[2] Uq = -ni = Uk = Uu = I; 1 ^ -k, k ^ 1; 

[3] Hq = -ni = -Uk = -Uu = -1; p -^ -k, k ^ p 

[4] Uq = -ni = -Uk = -n„ = 1; 1 - 



P 



[5] Uq = -ni ^ Uk ^ riu ^ -1; 
[6] Uq = ni = -Uk = Uu = -1; 



u, k^ -k, I 
u^ 1; 
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The delta functions and statistical factors have the same form for each of these six terms. Similarly, one can verify 
that the matrix element is the same in each case. The contribution to the collision term from p + q + I ^^ k + u 
processes is 
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where the matrix element squared of the diagrams shown in Fig. 1171 eauals 

\M\Il_^3[p + q + I ^ k + u] = ^ A-{p + q)[A-{k-l) + A-{k + u) + A-'{-l + u)] (124) 

+ A-{p + I) [A+ {q-k) + A+iq-u) + A-{k + u)] 
+ A- (p - k) [A+ (q + I) + A+iq-u) + A- (-[+ u)] 
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+ A-{p- u) [A+ {q + I) + A+{q-k) + A- {-1 + k)] 

+ A+{q + l)A-Ck + u) + A+(q- k)A-{~l + u) + A+{q- u)A-{-l + k) 
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FIG. 17: Scattering diagrams that contribute to the process p + g + /^^fc + uin the order they appear in Eq. H124t . 



The collision term for (^^ theory at the three-loop level, which represents the 2 <-^ 3 processes, is given by the sum 
of the contributions H120I123|) . Such a form of the collision term was postulated in j^, [ifl to reproduce the Kubo 
formula result of the viscosity coefficients within a linearized kinetic theory. 



XVI. FOUR-LOOP CONTRIBUTIONS 



At the four-loop level, we consider only 0^ theory. The relevant diagrams are the double sunset and two crossed 
versions of the double sunset, as shown in FigEl AH other contributions contain tadpoles and will not be discussed, 
as has been previously explained. 




FIG. 18: Four-loop contributions to the (^* self-energy. 



We follow the strategy formulated in Sec. lXIVBl and label internal lines of the diagrams in all possible ways that 
produce the same statistical factor. As previously, we define the operation 
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where we have used the definitions (|109|l and H116() . and additionally, we have introduced x = q+p + k + l + m. After 
summing up all central cuts of the diagrams shown in Fig 1181 we get the four- loop collision term 



e(po)v4(g)C4L-04(p) = —<^\A+{p + q + k) + A+{p + q + l) + A+{p + q + m) 
+ A'^(p + k + l) + A+(p + k + m) + A+(p + ^ + m) 



(126) 



+ A'^{k + l + m) + A+{q + k + l) + A+{q + l + m) + A+{q + k + m)\ >4L-,^* ■ 

We rewrite this result by introducing the delta function d^'^' {q+p+k + l-\-m — x) and an integral over x. As 
we have done previously, we will use Eq. ij9()(l and consider only positive energy contributions to the collision 
term. The sum over the remaining n's contains 32 terms. Five of these terms can be combined to give the process 
p + g^^fc + Z + m + x, ten give the process p-^q-\-l-\-m<-^ k -\- x, and ten give the process p + q + l'^m + k + x. 
The remaining seven terms are kinematically forbidden. Below, we give the collision terms representing these processes. 

The process p + q^^ k + l^m + x 



,+) , , , , 1 f d^q f dH f d^k f d^m f 

C,,^,4p + ,^l + k + m + x] = -J ^-^J ^^-^J j^-^J j^-^^j 



d^3 



{2TTfEkJ {2TTfE„J {2nfE, 
X {2'Kf6'^^\p + q-l-k-m-x) \M\li^_^i[p + q ^ I + k + m + x] 

X ^(1 + /(p)) (1 + /(q)) /(I) /(k) /(m) /(x) 

- /(P) /(q) (1 + /(I)) (1 + /(k)) (1 + /(m)) (1 + /(x)) 
where the matrix element squared of the diagrams shown in Fig. 1191 equals 



(127) 



\M\lL..[p + q^l + k- 



m -\- x\ = 
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A^{k + l + fh) + A+(g - fc - /) + A+(g - k - m) 
+ A+(g- [-m) + A+(p- k-l) + A+(p- k - in) 
+ A+{p - I- m) + A+(p + q-k) + A+{p + q-l) + A+{p + q-m) 



(128) 
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FIG. 19: Scattering diagrams that contribute to the process p + q*^k + l + m + x in the order they appear in Eq. 112811 . 



The process p + q + l + m^k + x 
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(+) z, 1 _ 1 [ Sq^ f dH f dhri f (fk f Sx 

C^^_^.[p + q + l + m^k + x\ - —J J^^^^J J^^^J j2^^fEZj WWJ WT^. 

X {2TTy6'-^'^{p + q + l + rh-k-S;) \M\l^_^4[p + q + I + m ^ k + x] 
X '{1 + /(p)) (1 + /(q)) (1 + /(I)) (1 + /(m)) /(k) /(x) 

- /(P) /(q) /(I)) /(m)) (1 + /(k)) (1 + /(x))^ , 
where the matrix element squared of the diagrams shown in Fig. 1201 equals 

A+(l + Th~k) + A+{q-k + l) + A+{q-k + m) 

+ A+(q + l + rh) + A+(p -k + l) + A+(p - fc + m) 

+ A+{p + 1 + m) + A+(p + q - fc) + A+{p + q + l) + A+{p + q + in) 



\M\lL_^4p + q + l + m^k + x] - ^ 



(129) 
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FIG. 20: Scattering diagrams that contribute to the process p + q + l + m^^k + x in the order they appear in Eq. IjlHOjI . 
The process p + q + l^m + k + x 



C 



(+) 

4L-. 



i[p + q + I ^^ m + k + x] = 



1 



dh 



dH 



d^r 



d^k 



dh 



2\i\J {2^YEqJ {2nYEiJ {2T:fE^ J {2nfEk J {2^YE^ 
X {2t: f 6^'^\p + q + I - rii -k - x) \M\li^_^i[P + 1 + ^ ^m + k + : 
X '{1 + /(p)) (1 + /(q)) (1 + /(I)) /(m) /(k) /(x) 

- /(P) /(q) /(I) (1 + /(m)) (1 + /(k)) (1 + /(x))] , 
where the matrix element squared of the diagranrs shown in Fig. 1211 equals 



(131) 



\'^\Il-4>4p + Q + >• ^m + k + : 



9_ 
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A+(/-fc-m) + A+{q-k + l) + A+{q-k-m) 



(132) 



+ A+{q + l- rh) + A+(p - fc + f) + A+(p - fc - m) 

+ A+(p + [- m) + A+(p + q~k) + A+{p + q + l) + A+(p + g - m) 
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FIG. 21: Scattering diagrams that contribute to the process p + q + l^^m + k + x in the order they appear in Eq. 113211 . 



The collision term of the 0^ theory at four-loop level, which represents the 2 <-^ 4 and 3 ^-+ 3 processes, is given by 
the sum of the contributions H127I129I131|I . As in the case of (j)'^ interaction, a form of the collision term beyond binary 
approximation was postulated in 9, 10] to reproduce the Kubo formula result of the viscosity coefficients within a 
linearized kinetic theory. While the contributions representing the 2 ^^ 4 processes coincide, the term corresponding 
to the 3 ^^ 3 interactions is absent in the effective transport equation 0, [T^l , as it is supposed to give negligible 
contribution to the viscosity [IJ. As already mentioned, the scattering 3 <-^ 3 was discussed in 8]. 

In performing the perturbative calculations of the collision terms in Sees. IXIIIIXVll we have used non-interacting 
Green's functions. However, the transport equation is of interest precisely because its solution is the full interacting 
distribution function. The terms on the left hand side of the transport equation clearly depend on this interacting 
distribution function. On the right hand side, we use the perturbatively calculated collision term and replace the non- 
interacting distribution function by the interacting one. We drop gradient terms in the expansion of the interacting 
distribution function about the non-interacting distribution, when these gradient terms appear within the collision 
term. The justification for this procedure is as follows: we consider the gradient and coupling constant approximations 
as independent, and choose to investigate higher order terms in the coupling constant approximation, while working 
at lowest order in the gradient approximation. The contributions we have calculated correspond physically to multi- 
particle production processes involving zero width quasi-particles. We remind the reader that finite width quasi- 
particles will give rise to contributions to the collision term that are kinematically forbidden for zero-width quasi- 
particles. For (j)"^ theory an example is a 1 <-> 3 process. In this paper we assume infinitely narrow quasi-particles (Eq. 
(|29|l '). and in this limit the effects of these processes will be negligible. 



XVII. SUMMARY AND OUTLOOK 



Kinetic equations usually take into account the interaction of particles with a mean field, and inter-particle pair 
collisions. However, when the system of interest is very dense, multi-particle interactions are expected to play a 
significant role in the system's dynamics. In addition, if a characteristic particle's kinetic energy is comparable to 
the particle's mass, particle production processes become important. Both multi-particle interactions and production 
processes occur in systems of relativistic quantum fields, such a quark-gluon plasma, when the energy density is 
sufficiently high. 

In this study, we have given the first systematic derivation of the transport equation of relativistic quantum fields 
which takes into account multi-particle and production process. Using the Schwinger-Keldysh approach, we have 
discussed scalar fields with cubic and quartic interactions. Mean-field phenomena are controlled by the one-loop 
tadpole contributions to the self-energy, and binary collisions correspond to the two-loop graphs. Multi-particle 
process occur at the three-loop level in cj)^ theory and at the four-loop level in 0* theory. Analysis of three- and 
four-loop diagrams is a technically very complex problem which requires special computational techniques. We have 
used the Keldysh representation because it provides a clean separation between propagators that correspond to real 
and virtual particles. In addition, we have made use of the cut structure of the collision term. Our calculations have 
been performed with the help of a program that uses MATHEMATICA symbolic manipulation software. 

Throughout this paper, we have assumed that quasi-particles are of zero-width and that their four-momenta are 
on the mass-shell. This assumption is crucial since some processes which occur at one- and two-loop level, like 1 <-> 2 
and 1 <-> 3 processes, are kinematically forbidden for on mass-shell particles. We have derived the explicit form of the 
contributions to the collision term of the transport equation corresponding to the processes: 2 ^-> 3 for the (j)^ model. 
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and 2 ^-> 4 and 3 ^^ 3 for the 4>'^ model. 

In this paragraph we clearly state the domain of applicability of our final results. Our derivation relies on several 
assumptions and approximations which are discussed in Sec. First of all, the system has to be homogeneous at a 
scale which is smaller than or comparable to the characteristic inverse momentum and the inverse effective mass, as 
required by the conditions H26I27I) . The system has to be weakly interacting so that the loop expansion is a legitimate 
approximation. This requires smallness of the coupling constants. Since we consider quasiparticles which are on the 
mass shell, the conditions (|27|l and H29|) have to be fulfilled. The latter is easily satisfied if the bare mass is much 
larger than the medium correction to the mass. A system of weakly interacting massive scalar fields which is close to 
equilibrium satisfies all of these conditions. 

In the future we plan to study further the multi-particle collision terms derived here. Although the integrals seem 
to be regular for massive particles, they require very careful analysis, as the transport theory of dense gases or liquids 
is known to suffer serious problems at the level of multi-particle interactions, see e.g. |42l |. The integral represnting 
the process 3 <-> 3 is of particular concern here as one particle can experience zero momentum transfer producing a 
singularity of the respective propagator. 

We intend to extend the work of Refs. |3, [ly| on the role of multi-particle processes in transport phenomena. 
However, we are going to start with the transport equations, which are derived here, and not with the Kubo formulas. 
We also plan to study how the system reaches chemical equilibrium due to particle production processes. The first 
step in this direction is to compute the interaction rate 



, ^ I / (IT) 

T[n ^ m] = - I C^+^[p + pi +P2 H Vn-i ^ gi + 92 H h Qm] , (133) 



dcf l_ f (Pp 

nj {2nfE, 

where C^^'[- ■ ■] is the gain contribution to the collision term corresponding to the process p + pi + P2 + ■ ■ ■Pn-i — * 
Qi + Q2 + ■ ■ ■ + Qm- Relative rates of the form r[2 — > 3]/r[2 -^ 2] or r[3 — ^ 3]/r[2 -^ 2] will give a measure of the 
importance of multi-particle processes. 
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